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Abstract. A library of 183 588 synthetic spectra based on Kurucz's ATLAS9 models is presented for the far red spectral 
interval (7653 - 8747 A). It is characterized by 3 500 K < T cff < 47 500 K, 0.0 < logg < 5.0, -3.0 < [M/H] < +0.5, 
< V r ot < 500 km s , £ = 2 km s _1 . The whole grid of spectra is calculated for resolving powers 8 500, 11 500 and 20000. 
A section of the grid is also computed for [a/Fe\=+0A and for microturbulent velocities and 4 km s . The library covers 
the wavelength ranges and resolutions of the two ambitious spectroscopic surveys by the ground experiment RAVE and the 
space mission GAIA. Cross-sections across the multi-dimensional data-cube are used to illustrate the behaviour of the strongest 
spectral lines. Interpretation of real data will have to include interpolation to grid substeps. We present a simple estimate of the 
accuracy of such a procedure. 
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1. Introduction 

Observations in the far red interval (7650-8750 A) are becom- 
ing increasingly used for determination of basic stellar param- 
eters. The advantages of this interval are its low sensitivity to 
reddening, high photon budget for late-type stars, good sensi- 
tivity of modern CCD detectors, and above all richness of its 
spectral features. 

The far red interval includes several important spectral 
lines (Munari 1999): (a) the K I infrared resonant doublet 
(7664.907 A, 7698.979 A), (b) the Na I non-resonant dou- 
blet (8183.256 A, 8194.821 A), (c) the O I triplet (8446.247, 
8446.359, 8446.758 A), (d) the Fe I multiplets 60, 401 and oth- 
ers with excitation potentials between 2.17 and 5.00 eV, (e ) the 
head of the Paschen series of hydrogen, and, most importantly, 
(/ ) the Ca II non-resonant infrared triplet (8498.018, 8542,089, 
8662. 140 A). The latter is present in all spectral types between 
B8 and M. A further advantage is that the interval of 8400 to 
8750 A is nearly free from telluric absorptions (Munari 1999). 

The Ca II triplet is a powerful diagnostic tool: the strength 
of Ca II lines depends on metallicity, but it is essentially in- 
sensitive to age of stellar population (Garcia- Vargas, Molla & 
Bressan 1998; Schiavon, Barbuy & Bruzual 2000; Vazdekis 
et al. 2003). These results have been used to study behaviour 
of Ca II in composite stellar populations and to attempt to 
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disentangle their age and metallicity. Composite systems in- 
clude normal elliptical galaxies (Molla & Garcia- Vargas 2000; 
Saglia et al. 2002; Cenarro et al. 2003), dwarf ellipticals 
(Michielsen, De Rijcke & Dejonghe 2003) and active galax- 
ies (van Groningen 1993; Diaz, Terlevich & Terlevich 1989; 
Nelson & Whittle 1999; Schinnerer, Eckart & Tacconi 2001; 
Marquez et al. 2003). Cenarro et al. (2002) studied the be- 
haviour of Ca II line-strength indices in terms of effective tem- 
perature, surface gravity, and metallicity. Their analysis was 
based on a set of 706 real spectra. A much larger synthetic spec- 
tral library presented in this paper can be used to complement 
these results on Ca II and other lines in the red domain. 

Two of the largest forthcoming spectral surveys are cen- 
tered on the far red spectral interval. GAIA is the cornerstone 
6 mission of ESA, approved for a launch around 2010. It is 
aimed at providing micro-arsec astrometry and ~ 10-band pho- 
tometry for ~ 10 9 stars brighter than V = 20. Brighter targets 
(V < 17.5) will also be observed by an on-board spectroscopic 
instrument operating in the 8480-8747 A wavelength range at 
a resolving power R = A/AA = 1 1 500 (Katz 2003). The other 
survey, Radial Velocity Experiment (RAVE), has just started in 
April (Steinmetz 2003). This is an international collaboration 
which uses the UK-Schmidt telescope at the Anglo Australian 
Observatory (AAO) equipped with a fiber-optic spectrograph to 
obtain spectra in a similar wavelength domain (8410-8750 A) 
and at a resolving power R ~ 8 500. The goal is to observe 
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Fig. 1. Graph of calculated spectra without alpha enhancement and with a micro-turbulent velocity of 2 km s~'. Metallicity is coded with 
different symbols which are plotted with small horizontal offsets for clarity. Each spectrum is calculated for 1 1 (T e g < 7000 K) or 14 (T tS > 
7000 K) different rotational velocities and at three different resolving powers (see text). 



35 million stars brighter than V = 16 at declinations suitable 
for AAO. The primary motivation for both surveys is to obtain 
stellar radial velocities to be used in studies of Galactic kine- 
matics. But a cross-correlation with a library of stellar spectra 



with known values of physical parameters can yield much more 
than radial velocity: in fact the effective temperature, gravity, 
metallicity (with certain element abundances), rotational veloc- 
ity, and the presence of different kinds of peculiarities can all be 
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determined. Therefore any software development for the GAIA 
mission as well as analysis of the first RAVE data needs to use 
an extensive library of stellar spectra in the far red domain. The 
same holds for individual spectroscopic studies in this wave- 
length range which are becoming ever more numerous owing 
to the preparation of the GAIA mission (see e.g. Munari et al. 
2001, Zwitter et al. 2003, Marrese et al. 2004). 



Table 1. Properties of models at different metallicities. The first three 
columns give values of metallicity, ff-enhancement coefficient and 
microturbulent velocity (£). The last two columns state if the "no- 
overshooting" approximation to the convection treatment and if new 
opacity distribution functions were used. 
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Fig. 2. As Fig. 1, but for models with alpha-enhancement ([a/Fe] = 
+0.4). All models were calculated for the metallicities [M/H] = -0.5, 
-1.0,-1.5, and -2.0. 
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Fig. 3. As Fig. 1, but for models with micro- turbulent velocities of 
and 4 km s' 1 . All models have solar abundances. 



Systematic observations of MKK standard stars in the far 
red spectral domain were presented by Munari & Tomasella 
(1999) and by Marrese, Boschi & Munari (2003). It is dif- 
ficult however to obtain high-signal-to-noise-ratio observa- 
tions to map all relevant combinations of stellar parameters 



in a uniform manner. So a parallel effort was launched to cal- 
culate a grid of synthetic spectra using ATLAS9 models from 
Kurucz. The first two papers (Munari & Castelli 2000, Castelli 
& Munari 2001) explored the grid in the temperature-gravity- 
metallicity space for non-rotating stars, assuming a micro- 
turbulent velocity of 2 km s _1 . Altogether 952 spectra were 
presented. 

These early works were aimed to coarsely but rapidly ex- 
plore an essentially unknown wavelength range to assist early 
planning and instrument design for GAIA. This phase is now 
over, and the community requirements are now moving toward 
data and analysis reduction pipelines. A much more extended 
and complete grid is now essential, and to provide one is the 
aim of this paper. The grid made available with this paper 
should meet the community requirements in testing algorithm 
coding for some time to come. The one to be used to ana- 
lyze the actual GAIA data when they will be finally assembled 
around 2015-2018 will be computed only in the next decade 
taking full advantages of the continuous advacements in the in- 
put physics, coding and atomic constants that will be reached 
by that time. The present grid is also timely presented to as- 
sist with analysis of RAVE spectra that are already routinely 
obtained at A AO. 

The present grid extends earlier calculations by adding 
more spectra, inclusion of stellar rotation and presentation 
of results at different spectral resolving powers. Furthermore 
a limited number of spectra corresponding to enhanced a- 
element abundances and different values of micro-turbulent 
velocity are presented. The database consists of more than 
183 500 spectra and is freely available in electronic form via 
CDS as well as via the ESA webserver. 
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2. Computation of synthetic spectra 

We computed synthetic spectra for the 7650-8750 A inter- 
val for almost all ATLAS9 grids of model atmospheres with 
different metallicities and micro-turbulent velocities available 
at the Kurucz web-site (http://kurucz.harvard.edu). We gener- 
ally used the NOVER models which differ from the previous 
Kurucz (1993a) models for the convection in that they were 
recomputed by Castelli with the overshooting option for the 
convection switched off (Castelli et al. 1997). When NOVER 
grids were not available, we used the Kurucz (1993a) mod- 
els. Because new ATLAS9 models based on updated Opacity 
Distribution Functions (ODFNEW models) have been recently 
computed (Castelli & Kurucz 2003), we used them for a few 
metallicities with T e g > 5000 K and for all the spectra with 
T eS < 4750 K. In fact, the main difference between the NOVER 
models and the ODFNEW models is the use of new Opacity 
Distribution Functions (ODFs) computed with TiO lines from 
Schwenke (1998) instead of from Kurucz, and including H2O 
lines, not considered at all in the previous ODFs. Also the 
ODFNEW models computed up to now are available at the 
Kurucz web-site (ODFNEW grids). Details of the adopted 
models are given in Table 1 . 

The synthetic spectra were computed with the SYNTHE 
code of Kurucz (1993b) at a resolving power of 500000. 
The source of atomic data was Kurucz & Bell (1995), the 
source of molecular data, except TiO, was Kurucz (1993c), 
while the source of TiO data was Kurucz (1999a), who sup- 
plied Schwenke's (1998) computations in SYNTHE format. 
Di-atomic molecules used in model computation are discussed 
in Kurucz (1992) and Kurucz (1993d). The source for H 2 
data is Partridge & Schwenke (1997), as distributed by Kurucz 
(1999b). The solar abundances are from Anders & Grevesse 
(1989), except for the spectra based on the ODFNEW mod- 
els. In this case the solar abundances are from Grevesse & 
Sauval (1998). Each spectrum computed for a given model at- 
mosphere was then broadened for several values of rotational 
velocity V rot and for three resolving powers: R =8 500, 1 1 500, 
and 20 000 by assuming a Gaussian instrumental profile. These 
resolutions were chosen because they correspond to the base- 
lined values for the RAVE survey, the GAIA mission and a 
typical Cassegrain-fed Echelle spectrograph, respectively. The 
resulting final synthetic spectra were resampled to 2 pixels per 
resolution element (R = 8 500, 1 1 500) or 2.5 pixels per resolu- 
tion element (R = 20000). 

3. Grid of synthetic spectra 

Ranges and steps for all seven basic parameters of the grid 
of synthetic spectra are given in Table 2. We adopt a com- 
mon convention of quoting metallicity and enhancement of a- 
elements in logarithmic units with respect to the solar values. 
The gravity is in logarithmic cgs units. Details of all calculated 
parameter combinations are given in Figures 1-3. Spectra are 
placed in gravity-temperature planes, with metallicity coded 
by a symbol type. Figure 1 covers the most numerous spectra, 
i.e. the ones with no a-enhancement and with microturbulent 
velocity of 2 km s _1 . The computed spectra cover the whole 



gravity-temperature plane except for hot low-gravity mod- 
els which are not radiatively stable. Low-temperature spectra 
(T e ff < 5000 K) were computed for a sparser set of metallici- 
ties due to large requirements of computing time. These spectra 
will be added online when completed. Fig. 2 corresponds to a- 
enhanced cases and Fig. 3 to those with a different value of 
microturbulent velocity. Note that each of the symbols actually 
corresponds to 1 1 (r eff < 7000 K) or 14 (r eff > 7000 K) spec- 
tra with different values of rotational velocity (see Table 1) and 
to three different resolving powers. 

All spectra are available as ascii files grouped into different 
directories according to their resolving power and temperature. 
The filenames are in a standard format identified in Table 3. So 
F765875v01 0r20000m05t05250g45k2nover.asc corresponds to 
a flux calibrated spectrum between 7650 and 8750 A, with 
V rot = 10 km s -1 , A/AA = 20000, [M/H] = -0.5, r eff = 
5250 K, log g = 4.5, £, = 2 km s~\ and no a-enhancement. 

The calculated grid is by far too large to present all of 
its properties here, so we explore only sample cross-sections 
across the grid. Figure 4 is a greyscale presentation of the 
spectra which were normalized to enhance line visibility. Each 
panel shows variation along one parameter axis, starting from 
a spectrum of a non-rotating K0 V type star. Note that all spec- 
tra were calculated in a wider wavelength domain, but only the 
8400-8750 A range is plotted for clarity. 

The temperature panel of Fig. 4 clearly shows the impor- 
tance of sharp Ca II lines for any radial velocity study. The 
panel is a textbook example of the expected behaviour of the 
Paschen lines and metallic lines. The metallicity panel illus- 
trates that the Ca II lines remain strong even at the lowest 
metallicities and the gravity panel shows their presence in all 
luminosity classes. The rotational velocity and resolving power 
panels show how the lines get smeared at high rotational veloc- 
ities or if observing at low resolving powers. 

The steps in the calculated grid are relatively small, but 
the coverage is not continuous. As an example, the step in 
temperature is 250 K (for T e s < 10000 K). This is larger 
than the baselined accuracy of temperature determination for 
both GAIA and RAVE surveys. So the grid will have to be in- 
terpolated to smaller steps. Figure 5 illustrates the errors in- 
troduced by a simple linear interpolation. At a certain grid 
point i with the parameter values /?, we compare the true 
synthetic spectrum S (pi) with the spectrum S ' obtained from 
a linear combination of spectra at neighbouring grid points: 
S' = fi-iS(pi-i) + f i+ iS(p M ). The weights and f i+l are 
optimized so that f[S(pi) - S'] 2 dA is minimal. The difference 
between the interpolated values of parameters p' and the true 
ones pi can then be expressed in units of a grid step: 

A = P ' - Pi - -/i+lQ^'+l ~ Pi) + fi-liPi-l ~ Pi) ^ 
Pi+1 ~ Pi-l (fi+l + fi-\)(Pi+\ ~ Pi-l) 

Figure 5 shows that linear interpolation is accurate to * 10 % 
of the grid step. Note that this is the worst case scenario, cor- 
responding to a reconstruction of the spectrum at the middle of 
the grid interval. Linear interpolation would be more accurate 
for spectra lying closer to one of the grid points. The results 
could be improved further by employing non-linear interpola- 
tion schemes. One may conclude that linear interpolation itself 
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Fig. 4. Cross-sections through the multidimensional [Temperature, Metallicity, Gravity (cm s~ 2 ), Rotation, Resolving power] data-cube of 
computed spectra. Parameter values correspond to a non-rotating K0 V type star (r eff = 5250 K, [M/H] = 0.0, logg = 4.5, V rot = km s -1 , 
£ = 2 km s -1 , A/AA = 20000), except for the parameter which is allowed to vary on a particular plot. All spectra were normalized by a single 
cubic spline fit to their upper envelope. The number of greyscale hues (identified at the top of the figure) is intentionally kept small, so as to 
allow reading the depth of individual spectral lines from the graph. Ions giving rise to the most prominent spectral lines are identified at the 
top. Synthetic spectra span a wider wavelength range and include variation of additional parameters (c-enhancement, different micro-turbulent 
velocities) not shown here for clarity. 

does not introduce errors exceeding 25 K in temperature (for ratio complicates their determination (Bailer-Jones 2003, see 
T e ff < 10 000 K), 0.05 dex in [M/H] or logg and 1 km s _1 in also Fig. 1 in Zwitter 2002). Also, spectra of real stars do not 
V ro t. Note that other errors are more important: degeneracy of correspond exactly to the synthetic spectra due to their pecu- 
parameter values fitting spectra with a limited signal to noise liarities (e.g. emission lines, varied abundances of individual 
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Table 2. General ranges of parameters of the spectra. Not all combinations were calculated (for details see Figs. 1-3). 
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Fig. 5. Errors introduced by linear interpolation of the grid. Each 
graph explores variation of one parameter, with the values of other 
parameters being fixed at: T eff = 5250 K, [M/H]=0.0, log£ = 4.5, 
V rot =0 km s" 1 , R = 20000. The ordinate is the difference between the 
value obtained from linear grid interpolation of the given parameter 
and the true value, expressed in units of the grid step. Note that linear 
interpolation yields values accurate to 1 10% of the grid step. 



Table 3. Filenames of individual spectra. Meaning of corresponding 
characters is given. 



character 


meaning 


1 


f: fluxed spectrum (erg s~' cm" 2 A steradian" 1 ) 




n: normalized to the continuum flux 


2-4 


starting wavelength (nm) 


5-7 


ending wavelength (nm) 


9-11 


rotational velocity V rot (km s" 1 ) 


13-17 


resolving power R 


18 


p: [M/H] > 0.0 




m: [M/H] < 0.0 


19-20 


10 x ABS([M/H]) 


21 


t: no a-enhancement 




a: [a/Fe] = 0.4 


22-26 


effective temperature r cff (K) 


28-29 


lOxlogg 


31 


micro-turbulent velocity £ (km s -1 ) 



4. Discussion and Conclusions 

The grid of synthetic spectra in the 7653-8747A interval pre- 
sented here is the most extensive so far. Altogether the grid 
contains 183 588 ascii files. They span three spectral resolu- 
tions corresponding to a typical Echelle spectrograph and the 
spectrographs of the RAVE and GAIA surveys. This should 
facilitate comparison of the results obtained with different in- 
struments. Spectra with other resolving powers (R < 20000) 
can be easily computed from the grid. 

The synthetic spectra we have computed can be used as 
templates for the determination of radial velocity, the primary 
goal of the RAVE and GAIA surveys. They also permit the user 
to derive the primary parameters of stellar atmospheres: tem- 
perature, metallicity, gravity and rotational velocity. The grid 
includes subsets of spectra computed for +0.4 enhanced abun- 
dances of the a elements and for different values of the mi- 
croturbulent velocity. It was shown that the grid can be easily 
interpolated to approximate spectra with parameter values be- 
tween the steps in the grid. The errors of parameters computed 
with a simple linear interpolation do not exceed 10% of the 
grid step. The computed spectra depend on the adopted values 
of individual element abundances. In particular the solar spec- 
trum computed with Anders & Grevesse's (1989) solar abun- 
dances for all the elements (Fig. 6) features too strong O I and 
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Fig. 6. Comparison of computed and observed spectra: (a) observed solar flux atlas resampled to R = 20 000 (Kurucz et al. 1984); (b) computed 
spectrum (T eB = 5750 K, [M/H]=0.0, logg = 4.5, V rot =2 km s"\ £ = 2 km s"\ R = 20 000) (c) same as (b) but for [M/H]=-0.5. All spectra 
are plotted in normalized flux. The topmost two spectra have been vertically offset for clarity. 



Fe I lines. This is consistent with a recently suggested lower 
oxygen solar abundance log(N(0)/N(H)) = -3.31 dex (Prieto, 
Lambert & Asplund 2001) and with a lower Fe I abundance 
log(N(Fe)IN{H)) = -4.5 dex (Bellot Rubio & Borrero 2002) 
as obtained from two and three dimensional hydrodynamical 
model atmospheres. This complements classical methods of di- 
agnostic line ratios for the far red spectral interval discussed by 
Munari (2002). 

Properties of Kurucz ATLAS model atmospheres in the 
context of the red spectral interval and the GAIA mission have 
already been reviewed by Nesvacil et al. (2003). Plez (2003) es- 
timated the potential of the MARCS models. Hauschildt et al. 
(2003) presented the computations of non-stationary and non- 
spherically symmetric atmospheres with significant amounts of 
dust. Finally Thevenin, Bijaoui & Katz (2003) reviewed the de- 
termination of chemical abundances from GAIA spectra. The 
common conclusion is that the observational potentials of the 
GAIA mission ask for improvements in the computation of 
synthetic stellar spectra. The RAVE collaboration is starting 
to obtain a very large number of high quality spectra in the 
far red spectral interval. They guarantee that this wavelength 
range is rapidly rising in importance. We are currently extend- 
ing the grid presented here to a larger wavelength interval. A 
continuous coverage from 250 nm to 1050 nm should facilitate 
a comparison of spectra obtained in other wavelength domains 
to huge observational sets coming from the RAVE and GAIA 
surveys. 
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